Conductance of DNA molecules: Effects of decoherence and bonding 
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The influence of decoherence and bonding on the linear conductance of single double-stranded 
DNA molecules is examined by fitting a phenomenological statistical model developed recently 
(EPJB 68, 237 (2009)) to experimental resuhs. The DNA molecule itself is described by a tight bind- 
ing ladder model with parameters obtained from published ab initio calculations (J.Am.Chem.Soc. 
127, 14894 (2005)). The good agreement with the experiments on sequence and length dependence 
gives a hint on the nature of conduction in DNA and at the same time provides a crucial test of the 
model. 

PACS numbers: 87.14.gk, 73.63.-b, 87.15.Pc 



I. INTRODUCTION 

Motivated by molecular electronics and by the pos- 
sibility to read out chemical and biological informa- 
tion electronically, transport through single double 
stranded DNA molecules is a focus of current research 
in nanoscience. The experimental situation is not par- 
ticularly clear. After the early investigations [T], in the 
last six years increasing consensus emerged that DNA 
double strands with 8 to 26 base pairs are conducting. 
However, this conclusion may be misleading, as it de- 
pends on the chemical potential of the contacts, on the 
temperature, and on the coupling of DNA to the envi- 
ronment. The experiments are performed either in water 
[2H1] or in the dry state ISj ■ Moreover, all experiments 
show sequence dependent conductance, which increases 
the more the (GC)-content in the sample molecule domi- 
nates over the (AT)-content. Sometimes ohmic behavior 
is reported, sometimes one finds an exponential decrease 
which is attributed to coherent tunneling. The change 
from coherent to ohmic behavior is usually attributed 
to the effect of decoherence caused by the water (e.g. 
0{S]), vibrational degrees of freedom [TU], or dynamical 
movement of the DNA bases [TTj . Most of the theoretical 
studies perform model calculations |12j , but some of them 
present hybrid methods combining ah initio and Molecu- 
lar Dynamics studies [B], sometimes also in combination 
with model calculations [TT] . 

The experiments clearly show that both the way, DNA 
is bonded, as well as the environment have an influence 
on the conductance. Here we present a calculation taking 
into account the effects of decoherence and bonding by 
extending a recently introduced phenomenological model 
[131 E]. It was originally developed for linear systems, 
and adapting it to the quasi-linear transport in DNA 
provides additional justification of the model. 



As DNA-Hamiltonian we use a realistic tight binding 
model (the extended ladder model) based on ab initio 
calculations [ISHTT] . 

We investigate several sequences, for which experimen- 
tal values for the conductance are known, in order to as- 
sess, how well they can be fitted by our model without 
changing the microscopic energy values, but just adapt- 
ing the Jour parameters fi and T describing the coupling 
to the electrodes, and p and rj describing the effect of 
decoherence (see Sections HI and IV I. 

The paper is organized as follows. In Sect. [IT] we 
present the extended ladder model describing the dou- 
ble stranded DNA. Sect. |IIl| is devoted to modelling the 
bonding to the electrodes according to the the experi- 
mental situation. In Sect. IIVI we extend the statistical 
model for decoherence such that it becomes applicable 
to DNA. In Sect. |V] we compare our results to the ex- 
periments 2-4 and finally, we give our conclusions in 
Sect. IVll 
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II. EXTENDED LADDER MODEL OF DNA 

The linear chain tight binding model |12) . where each 
site corresponds to a base pair, is widely used in the 
literature to describe DNA. In some works also the back- 
bone effects caused by the complementary strand and the 
sugar/phosphate mantle are taken into account (fishbone 
model) [IX. 

However, in order to account for arbitrary base se- 
quences considering the correct base pairing, the so called 
ladder model [T^l HH] is more appropriate. It consists 
of two coupled tight binding chains corresponding to the 
two strands, and each site represents a single base. There 
are calculations where backbone effects (here only due to 
the sugar/phosphate mantle) are considered in the ladder 
model, as well [T^ . 

Since, according to ab-initio |15H17| calculations the 
diagonal inter-strand transfer matrix elements can be 
more relevant than the intra-strand coupling, we will 
use the extended ladder model of DNA P^SHTT] sketched 




FIG. 1. The extended ladder model for DNA double strands 

[HE]- 



in Fig. [T] The corresponding Hamiltonian describing a 
double-stranded DNA that is N base pairs long is 

N 
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where c] „ creates a hole on strand m at the zth base pair 
with on-site energy e^^m. ti,mn and ii,i+i,mn are the base 
pair couplings and the hopping amplitudes between two 
bases of neighboring base pairs, respectively. The param- 
eters used in the calculations for the on-site energies and 
the base pair couplings are listed in Table |Tj According 
to |15j, the on-site energies calculated by density func- 
tional theory (DFT) depend on the flanking nucleobases, 
so that 16 different values were determined. For cq they 
vary between 7.890eV and 8.407eV, to give an idea. Here, 
we want to propose a universal model, which can predict 
conductances for any sequence of base pairs. Within the 
model class (fTl) a single value for the on-site energy of a 
nucleobase should be used. We simply take the average 
of the 16 values given in [15] for each nucleobase. 
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TABLE II. Hopping parameters (in eV) for the extended lad- 
der model, Fig.[l]p]. 
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FIG. 2. Density of states for holes, for 5'-(CG)^-3' 



TABLE I. The on-site energies and base pair couplings (in 
eV) . On-site energies are averages of the values calculated by 
DFT [15]. They correspond to the HOMO orbitals at the 
respective bases. 



following: The upper band at an energy — 8.2eV below 
the vacuum level and the lower one at — 9.7eV are both 
filled with electrons for an isolated neutral molecule. 



The hopping parameters (cf. Fig. [I]) used in the cal- 
culation are listed in Table |ll] We use the single-strand 
notation, listing only the sequence of a single strand (the 
other strand is determined due to the unique base pair- 
ing). Because of the directionality of the DNA strands, 
*5'-XY-3' ^ *3'-XY-5' = %'-YX-3' ^^^ ^ ^ Y. However, 
due to symmetry i5'_XY-5' ~ %'-YX-5' ^^'^ %'-XY-3' ^ 
^3'-YX-3' fo'^ all X,Y. 

Using the parameters given in Tables |l]and|n]the den- 
sity of states for holes of an infinite DNA double strand 
is shown in Fig. [2] We can see two bands, one for holes 
with an energy around 8.2eV and one around 9.7eV. Each 
band is split into two subbands of width sa 0.04eV. We 
prefer to describe the bands in terms of electrons in the 



III. BONDING MODEL 

Our aim is to examine the effect of the environment 
(decoherence) on the linear conductance of single double 
stranded DNA molecules. We focus on the experiments 
in [2H1] , where each gold electrode is coupled to a G-base 
on the 3'-end via a thiol group. We model this as shown 
in Fig. [3| assuming wide-band limit contacts attached 
only to the the 3'-ends of the strands. The self-energies 
due to the left (L) and right (R) contacts are 2N x 2A^ 
matrices, where all matrix elements are zero apart from 
the diagonal elements corresponding to (i,m) = (1,2) 
on the left and {i,m) = {N, 1) on the right. These ma- 





FIG. 3. Extended ladder model, where only the 3'-ends of the 
strands are attached to the gold contacts via thiol-linker like 
in the experiments of [2H1]- 



trix elements are equal to the imaginary number — ir/2. 
We will consider the case of weak coupling to the elec- 
trodes, by choosing F = 0.003eV. This is of the same 
order of magnitude as the weakest hopping parameters 
in the DNA molecule. 

The conductance of the DNA fragment depends also 
on the Fermi level of the contacts in relation to the DNA 
bands. Already a relative shift of O.OleV can significantly 
change the conductance. As the known values are hardly 
that reliable, we consider the chemical potential fj, of the 
contacts as a fit parameter in the following. It must be 
kept in mind that n denotes the chemical potential for 
holes. The one for electrons is given by /i^i = — /i. We 
concentrate on values of /id within the immediate neigh- 
borhood of the uppermost occupied band, i.e. around 
— 8.18eV. Hence we describe the bonding of the DNA 
fragment by two fit parameters, the self energy F and 
the chemical potential /t. 

In Sections |IV| and \V\ we will calculate, how the linear 
conductance of the experimentally investigated sequences 
depends on these two bonding parameters. They should 
be the same for all sequences in [2H1] , because the bond- 
ing was done the same way in all samples. We are going 
to use the recently developed phenomenological statisti- 
cal method [13 [H], which takes decoherence into account 
by two further fit parameters. 



IV. STATISTICAL MODEL FOR THE 
DECOHERENCE 

The model is based on a different physical picture com- 
pared to other phenomenological descriptions of decoher- 
ence (like the method of fictitious reservoirs (Biittiker 
probes) [19] or the method of [2D])- Whereas there de- 
coherence is continuously present in the sample, in our 
model it occurs only at stochastically distributed deco- 
herence regions, where phase information gets completely 
lost. 

We assume that decoherence can be described by lo- 
cal, stochastic events that couple a base pair (or part of 
it) to the environment. At those positions, phase coher- 
ence will be lost in both strands due to the interstrand 
couplings. In such a description the sequence consists of 
coherent sections separated by base pairs, on which deco- 
herence events take place. The positions of the decoher- 
ence base pairs are chosen at random with a probability 
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FIG. 4. The transmission function of a coherent double strand 
DNA molecule segment is calculated from the Hamiltionian 
given in Section jll] and the self-energies describing the con- 
nection to the decoherence base pairs at its ends (shaded). 



p giving rise to a particular decoherence configuration. 
The final results will be averaged over the different deco- 
herence configurations. The average decoherence length 
is a/p [13] , where a « 3.4A is the distance between base 
pairs, p is the first of the two parameters characterizing 
decoherence in the model. 

The coupling of a base pair to the environment im- 
plies that the local energy levels of the two bases will 
be broadened. In our model this broadening will be de- 
scribed by a parameter rj. As we furthermore assume that 
coherence is completely destroyed on these base pairs, it 
makes sense to attribute a local energy distribution func- 
tion to each decoherence base pair, which in general will 
differ from the Fermi distribution. Instead it will be a 
nonequilibrium distribution that emerges from the Fermi 
distributions fhiE) and /^.(-E) in the left, respectively 
the right contact via the transmission of charge carriers 
through the coherent sections. We expect 77 to be larger 
than the broadening ~ T/2 due to the electrodes and of 
the same order of magnitude as the width of the subbands 
inFig.[2](«0.04eV). 

The (n — l)st and the nth decoherence base pair can 
be regarded as "contacts" for the coherent section in be- 
tween. The electrical current (and hence the conduc- 
tance) can then be calculated according to the Landauer 
formula from the transmission function T'„_i „(£') of the 
coherent section and the energy distribution functions 
fn-i{E) and fn{E) at its boundaries. 



1^2- J dE r„_i^„(E) {fn-l{E) - fn{E)) . (2) 

The transmission functions for the coherent sections 
are calculated by applying the non-equilibrium Green's 
function method [21] (see Fig. [4]): 

Tn,7i+l(E) = 

4Tr[G„,„+i (E) Im S]„+i (£;)G1_„+i (£;) Im S„(E)] , (3) 
where 

Gn,n+l{E) = [E — Hn,n+1 — S„(J5) — E„+i(i5)] (4) 

is the Green's function of the segment between decoher- 
ence base pairs n and n + 1. It is a matrix of dimension 



2jn X 2j„, where j„ denotes the number of base pairs in 
the segment. The self-energy is a product of matrices 

(2j„ X 2i„) (2j„ X 2) (2 X 2) (2 x 2j„) 

(5) 
with the dimensions given in the second line. r„ is the 
coupling matrix connecting the segment to the decoher- 
ence base pair n (see Fig. B| . Its matrix elements can be 
identified from Fig. [l] The Green's function of the deco- 
herence base pair n contains the corresponding Hamilto- 
nian H„ and the coupling to the environment described 
in the simplest case by a constant imaginary self-energy 
— i?7 (times the 2x2 unit matrix): 



g^{E) = [E - H,, + iT^Y 



(6) 



As we assume complete loss of phase information at the 
decoherence regions, one can write down master equa- 
tions for the distribution functions, where the transfer 
rate of charge carriers between neighboring decoherence 
regions is proportional to the transmission function |13j . 
In the stationary state one obtains the following system 
of coupled linear equations: 

Tn-l,niE)[fn-l{E)-fniE)] = Tr,,n + l{E)[UE)~ f,,+ ^{E)]. 

(7) 
According to Eq. (pi), equations (l7| imply current con- 
servation. They can be solved analytically 13J, using 
the boundary conditions fo{E) — fhiE) and fM+i{E) — 
fR{E), where M denotes the number of decoherence base 
pairs. Inserting the solution into Eq. (pi) (e.g. for n = M) 
gives the current 



J '^'^ is ^"-+ J 



{fL{E)-Jn{E)). (8) 



This leads to the linear conductance 



dL 



cq 



dE 



(9) 



where Gq = 2^ is the conductance quantum (for the 
two spin channels), /eq is the Fermi function and /i is 
the equilibrium chemical potential of the contacts. For 
the results presented in the following, this expression has 
been averaged over typically 2000 independent decoher- 
ence configurations. 

It is worth noting that temperature enters the result 
via several parameters. First, there is the explicite tem- 
perature dependence of the Fermi function. Moreover, 
the chemical potential /i and the decoherence parame- 
ters p and 7] implicitely depend on temperature. 

Given the microscopic energy values in Tables |l] and 
ITT] and the temperature, the model has four parameters. 
Two of them, /j, and F, describe the coupling to the con- 
tacts. The other two are the decoherence parameters, p 
and rj. As they describe a coupling to the environment, 
they may differ for dry and wet DNA, and they may also 
depend on the type of base pair. 



V. RESULTS 

In this Section we present our results on the lin- 
ear conductance for the sequences investigated exper- 
imentally in [2Hlj. The experimental data consist of 
conductance values for 5'-(CG)m-3' samples [H E] with 
four different lengths (numbers of base pairs (uBF) are 
2m = 8,10,12,14), and for 5'-CGCG-(A)„(T)„-CGCG- 
3' samples [H [3] with three different lengths (uBPs are 
8-|-2n = 8, 10, 12), where n = corresponds to the short- 
est 5'-(CG)m=4-3' sample. 

We took the effect of finite temperature into account 
only in the Fermi functions of the electrodes. Since the 
experiments were performed at room temperature, we 
fixed ksT = 0.0255eV. It is interesting to note, that in 
[4] experimentally unobservable temperature dependence 
was reported in the temperature range from 5 to 65 Cel- 
sius. 

It is a hard task to investigate the whole four di- 
mensional parameter space of the model presented in 
Sect. |TVJ One needs to have reasonable guesses for the 
parameter values F, /x, p, and rj. It turns out that 
the chemical potential fi of the electrodes is the most 
crucial parameter. Here we present data for four val- 
ues, fi — 8.07eV (/Zoi slightly above the upper band), 
fi = 8.18eV (/ioi within the upper band), and fi — 8.25eV, 
respectively /x = 8.26eV (both /id slightly below the up- 
per band of the isolated molecule) . The coupling F to the 
electrodes is assumed to be of the same order of mag- 
nitude as the weakest hopping parameters in the DNA 
molecule. It is fixed to the value F = 0.003eV in all fig- 
ures apart from the last one. As explained in Sec |IV[ we 
regard values of r] slightly larger than 0.04eV (but still 
of the same order of magnitude) as reasonable. Finally, 
we guess that at room temperature a hole cannot travel 
ballistically further than perhaps two base pairs. Thus 
we choose p — 0.5. These guesses will be confirmed by 
the results. 

Fig. [5] shows the measured conductances of four 
molecules of the type 5'-(CG)m-3', plotted versus the 
inverse of their length L — Ira (number of base pairs 
(uBP)) [21 13]. Within the experimental error bars we 
obtained four excellent fits by keeping p — 0.5 and 
F = 0.003eV fixed and optimizing 77 for each value of 
the chemical potential. The statistical errors of the cal- 
culated conductances are not shown in this figure for the 
sake of clarity. They are always less than half of the 
experimental errors. In the cases, where the chemical 
potential lies outside the band, the value 77 — 0.05eV led 
to good fits for all four data points. However, if the chem- 
ical potential lies inside the band, a good fit requires a 
value of ry = leV, which seems unreasonably large. 

Plotting the conductance as a function of inverse 
length suggests ohmic behavior. This is misleading, 
though, because a linear extrapolation would lead to van- 
ishing conductance at a finite length. In Figs.[6][7| and [8] 
we plot the same data in a semilogarithmic way, showing 
that an exponential length dependence cannot be ruled 
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FIG. 5. Calculated conductances for 5'-(CG)m-3' versus the 
inverse of the number of base pairs (uBP) 2m for parameter 
combinations {^J. = 8.07eV, 77 = 0.05eV) (red), {/i = 8.18eV, 
77 = l.OeV) (blue), {/i = 8.25eV, 77 = O.OSeV) (green), and 
ifi = 8.26eV, 77 = 0.05eV) (orange). T = O.OOSeV, p = 0.5 
and fcflT — 0.0255eV are the same for all curves. The exper- 
imental values i2j i4j are shown by ■ together with the error 
bars. 



out. 

These Figures also show two more experimental data 
points, which belong to molecules 5'-CGCG-(A)„(T)„- 
CGCG-3' with n — I and n — 2. We calculated the cor- 
responding conductances using the same parameters as 
obtained in Fig. [5] for the other molecules. Only the on- 
site energies and hopping parameters were adapted ac- 
cording to the Tables H] and [n] For the cases fi = 8.07eV, 
/i = 8.25eV, respectively fj, — 8.26eV excellent fits for the 
two new data points were obtained (see Figs. [Gland l8|. 
However, for 11 — S.lSeV, the new data points cannot be 
fitted with the old parameters. This casts further doubt 
on the validity of this parameter set, in addition to 77 
being unreasonably large. 

In the following we examine, how sensitively the cal- 
culated conductance for the longest of the 5'-(CG)m-3' 
samples {m = 7) depends on the parameters 77, p, and 
r. The purpose is to show that the values leading to 
the fits in Fig. [5] do not offer much freedom of choice, 
when varied individually. The chemical potential is fixed 
at ^ = 8.07eV, /z = 8.18eV, respectively /i = 8.25eV, 
corresponding to the fit curves close to the upper end of 
the error bar for this sample. It should be kept in mind, 
that a change of parameters that moves this conductance 
closer to the experimental value, will in general lead to 
worse fits for the other three samples in Fig. [5] 

In Fig. |9] the calculated conductance is plotted as a 
function of 77 while keeping p = 0.5 and F — O.OOSeV 
fixed on the values, for which the fits in Fig. [5] were ob- 
tained. For fi = 8.07eV (red) and /^ — 8.25eV (green) 
the conductance does not depend sensitively on 77: One 
would get acceptable fits for values in a large interval be- 
tween O.OleV and 0.5eV. This is true for a molecule that 
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FIG. 6. Conductance data for 5'-(CG)™-3' (D) [1|1] and 5'- 
CGCG-(A)„(T)„-CGCG-3' (0) [lEl plotted semilogarithmi- 
cally versus the number of base pairs (nBP), 2m, respectively 
8 -I- 2n. For a value of /i = 8.07eV, our model fits both data 
sets, using the same parameters F = O.OOSeV, p = 0.5, and 
ri = 0.05eV. keT = 0.0255eV. 
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FIG. 7. Same experimental data as in Fig. [6] For a value 
of /i = 8.18eV, the parameter values fitting 5'-(CG)m-3', F = 
O.OOSeV, p = 0.5, 77 = l.OeV, do not fit the data for the chains 
with AT-inclusions. ksT = 0.0255eV. 



is 14 base pairs long. The optimization of the 77-value 
has to rely on the length dependence of the conductance. 
These data are crucial to narrow down the interval for 
this model parameter. By contrast, for /i = 8.I8eV the 
conductance depends very sensitively on 77, allowing ac- 
ceptable fits only in a narrow interval around leV and 
perhaps around a second value much below O.OleV. Both 
values are outside the range, which we regard as plausi- 
ble. 

In Fig. [To] the calculated conductance is plotted as a 
function of p while keeping 77 = O.OSeV respectively 77 — 
leV and F = O.OOSeV fixed on the values giving the fits in 
Fig. [5] For the chemical potentials outside the band (/i = 
8.07eV (red) and fj, = 8.25eV (green)) the sensitivity is 
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FIG. 8. Same experimental data as in Fig. [6] For values 
between /i = 8.25eV (green) and /i = 8.26eV (orange), our 
model fits both data sets using the same parameters V = 
O.OOSeV, p = 0.5, r] = 0.05eV. ksT = 0.0255eV. 




FIG. 9. The sensitivity of the calculated conductance on rj. 
The functions G{r]) are shown for /i = 8.07eV (red), p = 
8.18eV (blue), and ^ = 8.25eV (green) for 5'-(CG)7-3'. F = 
O.OOSeV and p — 0.5 are fixed. The experimental value is 
indicated by the horizontal full line, its errorbar by the dashed 
ones. fcsT = 0.0255eV as in the experiment. 



high: Only for p-values in the interval between 0.3 and 
0.5 one gets acceptable fits. However, for the chemical 
potential within the band (/i = 8.18eV) the sensitivity is 
low: Considering only the sample with 14 base pairs, any 
p-value between 0.3 and 1 would be allowed. 

In Fig. [TT] the calculated conductance is plotted as a 
function of F while keeping p — 0.5 and i] fixed on the 
values giving the fits in Fig. [5] For all three chemical 
potentials the conductance depends sensitively on F. Ac- 
ceptable fits are obtained between values somewhat be- 
low O.OOleV and O.OOSeV, justifying our weak coupling 
assumption. There may be a second fit interval with very 
large values above leV, which we discard as unplausible, 
because they would imply a coupling of the molecule to 




FIG. 10. The sensitivity of the calculated conductance on p. 
The functions G{p) are shown for jj, = 8.07eV, 77 — 0.05eV 
(red), II = 8.18eV, r? = l.OeV (blue), and /i = 8.25eV, 
77 = 0.05eV (green) for 5'-(CG)7-3'. F = O.OOSeV. The ex- 
perimental value is indicated by the horizontal full line, its 
errorbar by the dashed ones. ksT = 0.0255eV as in the ex- 
periment . 
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FIG. 11. The sensitivity of the calculated conductance on F. 
The functions G(r) are shown for p. = 8.07eV, r) = 0.05eV 
(red), /i = 8.18eV, r) = l.OeV (blue), and ^ = 8.25eV, t) = 
0.05eV (green) for 5'-(CG)7-3'. p = 0.5. The experimental 
value is indicated by the horizontal full line, its errorbar by 
the dashed ones. fcaT = 0.0255eV as in the experiment. 



the electrodes that is better than any coupling between 
the base pairs within the molecule. 



VI. CONCLUSIONS AND OUTLOOK 

We have proposed a simple model in order to calculate 
the conductance of DNA-double strand molecules for any 
sequence of base pairs. It is a tight binding model with 
on-site energies and charge transfer integrals taken from 
the DFT calculations in HIl. 



The on-site energies correspond to the energy dif- 
ference between the vacuum level and the highest oc- 
cupied molecular orbital (HOMO) of the nucleobase. 
In the 5'-(CG)oo-3' double strand they give rise to a 
lower band about 9.72eV, and an upper band about 
8.18eV below the vacuum level. Both of them arc filled, 
if the molecule is neutral. The energy to create an 
electron-hole pair in such a molecule is of the order of 
min(LUMOc,LUMOG) + 8.18eV w 3 - 4eV [J, where 
LUMOc/G denotes the energy of the lowest unoccupied 
molecular orbital of C respectively G. Actually the dif- 
ference between the band centers overestimates the band 
gap, but as the band width is of the order of 0.05eV, 
the conclusion is still correct that thermal activation of 
charge carriers at room temperature can be ruled out. 

The situation changes, if the DNA-molecule is brought 
into contact with gold electrodes. Let us first discuss the 
case, where the chemical potential of holes in the elec- 
trode is /i = 8.07eV (in other words, electronic states are 
filled up to the energy — 8.07eV). At room temperature 
there will be thermally activated holes in the electrode 
down to the HOMOq level. They can be injected into 
the molecule and carry a current through the molecule. 
In this sense the electrode has a similar effect as accep- 
tors in a p-doped semiconductor. This works only, if the 
molecule is short enough that no localized space charges 
form at the contacts. In summary, for /i = 8.07eV the 
molecule is a hole conductor, which remains neutral when 
contacted. 

If the chemical potential of holes in the electrode is 
H = 8.26eV, we obtained fits of the experimental data of 
similar quality. In this case, however, the Fermi energy of 
the gold electrodes lies below the uppermost filled band of 
the DNA molecule. Again we assume that the molecule 
is short enough that space charges extend all the way 
from one side to the other. Then the upper band will 
be depleted, the molecule will carry a positive charge of 
2N elementary charges, where N is the number of base 
pairs. In the electrodes there will be thermally activated 
electrons reaching up to the now empty HOMOq level. 



In this case the electrodes act like donors in an n-doped 
semiconductor. In summary, for fj, = 8.26eV the molecule 
would become charged when contacted, and it would be 
an electron conductor. 

The third parameter set, for which we could find an 
excellent fit of the conductance of 5'-(CG)„-3' double 
strands, had fi — 8.18eV, which lies within the upper 
occupied band of the DNA molecule. This will lead to 
a depletion of the band, which is only partial, so that 
no thermal activation of charge carriers is needed. In 
this sense the DNA molecule behaves like a metal, pro- 
vided it is shorter than the screening length. However, 
the fit requires a decoherence induced level broadening 
T] — leV, which is unphysically large. Moreover, the 
same parameter set is incompatible with the conductance 
data for the mixed 5'-CGCG-(A)„-(T)„-CGCG-3' double 
strands. Therefore, this scenario can be ruled out. 

Most of the experimental and the theoretical examina- 
tions refer to hole conduction. Therefore we can conclude 
that the first parameter set (/i = 8.07eV, 77 = O.OSeV, 
r = O.OOSeV, and p = 0.5) is the most plausible one to 
describe the DNA molecules considered. 

Summarizing, we extended our phenomenological 
model [131 [13] for describing the effect of decoherence 
in double stranded DNA molecules. We could fit all of 
the experimental data of [2H3] without changing the mi- 
croscopic energy values and by using the same parameter 
set of the model. For other experimentally investigated 
sequences [5J H] work is in progress. 
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